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Abstract 

By evaporating a drop of lipid dispersion we generate the myelin morphology often seen 
in dissolving surfactant powders. We explain these puzzling nonequilibrium structures 
using a geometric argument: The bilayer repeat spacing increases and thus the repulsion 
between bilayers decreases when a multilamellar disk is converted into a myelin without 
gain or loss of material and with number of bilayers unchanged. Sufficient reduction in 
bilayer repulsion can compensate for the cost in curvature energy, leading to a net stability 
of the myelin structure. A numerical estimate predicts the degree of dehydration required 
to favor myelin structures over flat lamellae. 

1 Introduction 

The phospholipid molecules that constitute the chief component of cellular plasma membranes 
spontaneously self-assemble to form bilayers when dissolved in water. These bilayers tend 
to stack to form multilayers, also known as multilamellae. At high lipid concentrations, 
the stable morphology consists of flat, stacked bilayers with quasi- long-range order ^|2HH]- 
Yet, under certain nonequilibrium conditions, the bilayers bend to form long-lived spherical 
multilamellae known as onions (3J|!)] or nested cylindrical tubes called myelin figures or simply, 
myelins (see Fig ^) 5, El El El, despite the additional penalty in curvature energy. The 
formation of myelins has been a mystery since their discovery more than 150 years ago[BlEl- 
Myelins offer a variety of potential applications such as encapsulation and controlled delivery 
of drugs 10]. Thus a better understanding of their formation and structure is desired. 

In this paper we identify a type of constrained equilibrium that leads naturally to myelin 
formation. The simple geometrical mechanism proposed here offers a plausible and general 
explanation for the longstanding puzzle of myelin formation and has explicit implications 
about their internal makeup. We note that in contrast to the Hij-to-L a -to-Hu reentrant 
phase transitionjllj that involves topological changes, the formation of myelins in our system 
involves only geometrical changes. 

This paper is organized as follows: In Section [21 we outline our experimental results of 
myelin formation. Based on the experimental observations, we propose a model in Sections 
121 and 21 to explain the formation of myelins in our experiment. In Section \5\ we investigate 
the effect of water permeation on our model. In Section we discuss the implications and 
limitations of our model. Finally Section [7| concludes our work. 

*email: jhuang2@uchicago.edu 



Figure 1: Formation of myelin figures in a typical contact experiment. A lump of concen- 
trated dilauroyl phosphatidylcholine (DLPC) (bottom) is contacted with water (top) at 27°C. 
Closely packed myelins form at the interface and grow toward the water-rich side. The ex- 
perimental approach has been described in detail in many references such as 00 EI- The 
bar at the lower left corner represents 100/im. 

2 Experiment 

The main goal of this paper is to present our model of myelin formation. Therefore we 
only summarize our experimental results that are pertinent to the model. More complete 
description of the experiment will be given elsewhere |12j. 

In our experiment we observe a drop of dilute suspension of dimyristoyl phosphatidyl- 
choline (DMPC) bilayer vesicles in water as it evaporates ^2] (FigEJ). The drop is heated to 
27— 30°C, above the chain-melting temperature of DMPC (~ 24°C)[T3]. Thus the bilayers are 
in the fluid state (L a phase) pQ. The evaporative flow creates a ring deposit of concentrated 
lipid around the drop's perimeter |14). From this deposit many pancake-like multilamellar 
disks develop and grow inward. As the disks grow, they undergo a shape transition to form 
myelins. Experimental evidence suggests that materials, i.e., lipid and water enter the disk- 
myelin complex mainly via its root embedded in the dry deposit region (FigEJ) IB]- If 
evaporation is halted, growth stops and the myelin is resorbed into its parent disk; but the 
myelin's cylindrical morphology is retained during the resorption. This suggests that the 
disk- myelin complex is in quasi-equilibrium^3], and hence free energy analysis is applicable 
to the disk-to-myelin shape transition. 

3 Model Definition 

In this section and the next one, we present a theory to account for the myelin formation, 
i.e., the disk-to-myelin shape transition, observed in our experiment (Fig[2J). We will show 
that the formation of myelins can be favorable because their cylindrical form results in a 
larger bilayer repeat spacing (i.e., separation between bilayers) and hence lower inter-bilayer 
repulsion than flat multilamellae. 

We postulate that (a) the number of bilayers is unchanged during the shape transition, 
and (b) the bilayers can exchange materials, i.e., lipid and water, freely to achieve quasi- 
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Figure 2: Optical micrograph of disk- myelin complex. Myelins emerge from a multilamel- 
lar disk at the perimeter of a water drop containing dilute dispersion of DMPC bilayer 
vesicles [T2], A: overview of the drop. B: myelin. C: multilamellar disk. D: initial contact 
line. E: dry DMPC deposit. F: current contact line. G: a band of aggregated vesicles. At 
the moment of the photo, only the long curved myelin was growing at about 0.5/iin/s. At 
later times, the long myelin drifted into place beside the straight myelins and another myelin 
pinched off from the tip of the disk, repeating the cycle. The bar at the bottom represents 
20/iin. 
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Figure 3: Coss-sectional view of a multilamellar disk (B) and a myelin (C). Both structures 
are composed of N uniformly spaced bilayers. Figure shows the special case of iV = 3. 
Oval inset (E) shows placement of lipid molecules (L) making up bilayers. The symbol S 
represents the mid-surface of a bilayer. The total bilayer area A is measured at S. In this 
work we neglect the thermal undulations of bilayers (see Section EJ) 23 , 24 . Therefore A can 
be determined simply from the geometry of the disk or myelin. The bilayer repeat spacing 
Dd or D m = d = d w + di, where d w is the water layer spacing and di the bilayer thickness. If 
the inter-bilayer pressure is sufficiently low, d\ is roughly constant and therefore the change 
in d is mostly due to that in ^(Fig.l of |20]). 



equilibrium^]. The number of bilayers, denoted by N, is unlikely to change on the time 
scale of the shape transition, which is about 1 second. The bilayers may exchange materials 
via defects such as screw dislocations !^! IT71 ITS] . 

The free energy includes a curvature energy favoring flat bilavers[l9] and a repulsion 
between bilavers |2*Ul 1211 122 j . The curvature energy of a bilayer of area A takes the form of 
(Fig© 

^ J A dA( Cl + c 2 -c f, (1) 

where k c is the bending stiffness, c\ and C2 are the principal curvatures, and cq is the spon- 
taneous curvature |19j. The bending stiffness 

k c = 1.2 x 1(T 19 J (2) 

for DMPC hi layers 2") . (However, Ref.[26; gives a smaller value of k c .) In our experiment we 
expect Co ^ 0. Postulate (a) implies that the topology of the disk-myelin complex does not 
change in the shape transition. Therefore we can neglect the Gaussian curvature contributions 
according to the Gauss-Bonnet theorem |27j. With (^Q) the curvature energy of a myelin 
consisting of N concentric, uniformly spaced bilayer cylinders is given by 

Kc-^-ln(^V (3) 



D m \ R, c 

where L is the myelin length and D m the bilayer repeat spacing (FigJSJ); R Q and Ri are the 
radii of the outermost and innermost cylinders respectively. The total bilayer area 

A = TrN{Ri + R )L. (4) 
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By construction 

R = Ri + (N-l)D m . (5) 
In this paper we set (see Section EJ) 

Ri = ^ (6) 

for convenience, and we only consider myelins of large aspect ratios, i.e., L ^> R Q , so that the 
end caps are negligible. 

The repulsion between bilayers is determined experimentally: The authors of [20] found 
that their pressure data were well represented as the sum of an exponentially falling hydra- 
tion force plus a van der Waals attraction (Table 1 of (SOI)- Our estimates of inter-bilayer 
interaction energy use this functional form. Specifically, given water layer thickness d w and 
bilayer repeat spacing d (FigEJ), the inter-bilayer pressure 

P = Pi+P 2 , (7) 

where the hydration pressure 

Pi{d w ) = P exp 

and the van der Waals attraction 



d w 

T 



(8) 



P 2 (d w ) = ~ 

07T 



1 2 i i 



(9) 



_dl d? {d + d w ) 3 _ 
The inter-bilayer interaction energy per unit bilayer area is therefore equal to 

d{d w )P{d w ). (10) 

o 

In this paper the bilayer thickness di is taken to be constant (see Fig|3J). The pressure 
Po = 10 9 ' 94 dyne/cm 2 and A = 0.26nm for DMPC bilayers, and the Hamaker coefficient H is 
set to 10 _20 J in this paper j2U| 128 ] . When the bilayers are curved, this energy is in principle 
altered. In Section we argue that such corrections are negligible in the case of interest. 

Postulate (b) proposed at the beginning of the section implies that the important external 
parameters for determining the free energies are the amount of lipid and that of water in the 
disk-myelin complex. Under moderate external pressure, the bilayer thickness di is approx- 
imately constant and, moreover, both lipid bilayers and water are virtually incompressible 
(FigEJ) 20 . This means we can describe the amount of lipid by the total area A of all the 
bilayers, which is measured at the bilayer's mid-surface (S in FigEJ). Similarly, the total 
amount of lipid plus water is equivalent to the total volume V of the complex. 

Given the total bilayer area A, volume V, and the number of bilayers N, the disk-myelin 
complex is internally confined. In this case the lipid concentration A/V determines the bilayer 
repeat spacing and thus the inter-bilayer repulsion, as shown in the next section. 



4 Instability of the Disk 

We now illustrate how a sufficiently large lipid concentration, i.e., area-to- volume ratio A/V, 
can make a multi-lamellar disk unstable to myelin formation. 

We consider a large but thin disk composed of N nested bilayer disks, with uniform bilayer 
repeat spacing Dd, volume V, and total bilayer area A. Because the disk is large and thin, 
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its rim part is negligible. Figure El shows the example with N = 3. If the area of each bilayer 
is denoted by A±, then the total area A = 2NA\ and the total volume V = (2N — 1)DjAi- 
Thus the repeat spacing 

D <=( 1 + 2ivb)T 

This Dd determines the inter-bilayer interaction energy (see (|1U|0. For these flat bilayers, the 
curvature energy © is arbitrarily small. 

To show the instability of this disk, we convert it into a myelin composed of N concentric 
and uniformly spaced cylindrical bilayers, keeping both V and A fixed (Fig©. Using 
© and © the bilayer area A = -KLN 2 D m while the myelin volume V = 7rR 2 L = irL(N — 
l/2) 2 D 2 n . Thus Vj A = (1 - l/2N) 2 D m so that the myelin repeat spacing 

1 \ 2 V 

D m = (1 + 



2N-1J A 

This D m determines the inter-bilayer interaction energy of the myelin (see (|l(Jjl ). With © 
and © its curvature energy © can be written as 

irL 

K c —ln(2N-l) (13) 

Because the myelin has a larger repeat spacing, it has a lower repulsion energy than the disk. 
If the decrease in the repulsion energy can compensate for the myelin's curvature energy 
(see (|1U|) and (|13j) ^. the myelin, rather than the disk, becomes the thermodynamically more 
stable morphology. Here we define a threshold inter-bilayer pressure Pth, above which the 
disk becomes unstable against myelin formation. Namely, Pth is the pressure at which the 
decrease in the inter-bilayer repulsion energy due to the disk-to-myelin transition is equal to 
the curvature energy of the myelin. 

Although the fractional change in spacing 

D "- D * = -J— „ J_ (14) 
D d 2N - 1 2N y ' 

is tiny when iV S> 1, the reduction in the total repulsion energy is proportional to N, and 
therefore the effect of spacing increase can still be significant, as demonstrated by the example 
below. 

Since the inter-bilayer pressure increases rapidly with decreased repeat spacing (J7J), a 
multilamellar disk with modest dehydration can easily have enough repulsion to become 
unstable. A typical disk observed in the experiment should be somewhat dehydrated because 
part of it is in the dry deposit region (FigJ2J. The disk thickness is observed to be about 2.5 
/um|12|. Given the equilibrium spacing 20 , we infer the number of nested bilayer disks to 
be approximately 200. 

Taking N = 200 and using the measured bending stiffness ((2j) and inter-bilayer repulsion 
(|7|) for these DMPC bilayers [2*0]. our model predicts that the energy of the disk becomes 
larger than that of the myelin when the dehydration exceeds 2.1%, i.e., a 2.1% decrease 
in the repeat spacing from the equilibrium value of 6.22nm[2l7]. The inter-bilayer pressure 
P = Pth ~ 0.32atm at such dehydration 0. The myelin converted from this disk has a repeat 
spacing about 1.9% less than the equilibrium value and an inter-bilayer pressure P of about 
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0.3atm, with the van der Waals attraction P2 — 0.013atm 0. Since the difference between 
the repeat spacings, Dd and D m , is tiny, the myelin diameter is nearly identical with the 
thickness of the parent disk; this is consistent with our experimental observations 12.. Here 
we assume that the threshold pressure Pth is sufficiently low so that the bilayer thickness di 
is unchanged during the shape transformation (FigEJ). With the experimental data shown 
in Fig.l of |20| . our result of Pth ~ 0.32atm suggests that a constant d\ is indeed a fair 
approximation . 

Based on our model we propose the following scenario to explain the myelin formation in 
our experiment (FigJ2J: As water evaporates, lipid and water enter the disk, and hence both 
A and V increase. Because the contact line gradually moves inward due to water evaporation, 
the disk becomes more and more dehydrated, i.e., the lipid concentration A/V increases. An 
increase in A/V leads to a decrease in the repeat spacing Dd (jllj) and thus an increase in 
the inter-bilayer repulsion (J7J). The repeat spacing keeps decreasing until a threshold, below 
which the disk is unstable and a myelin with a larger spacing grows out of it in order to lower 
the free energy. 



5 Effect of Water Permeation 

The above simple example has a large pressure difference ~ 0.3atm between the inside and 
outside of the disk-myelin complex. This pressure might induce sufficient water permeation 
through the bilayers to reduce the dehydration and restore the stability of the disk morphol- 
ogy, thus invalidating our model. In the following we perform a self-consistency check to show 
this is not the case: In practice the system should be in a steady state, where the pressure drop 
occurs across a significant fraction of the complex. This suggests that the pressure difference 
across a bilayer in the complex of Figure[2]is about latm/iV ~ latm/200 = 5 x 10~ 3 atm. The 
water permeability coefficients p w for typical phosphatidylcholine bilayers are known to be 
30-100 /im/s(29]. Assuming the validity of the solubility-diffusion mechanism |30j. the volume 
flux of water across a bilayer of area A2 is given by 

j - v Y^ A p 

where AP is the pressure difference across the bilayer, V w = 18cm 3 /mole, the partial molar 
volume of water, R the gas constant and T the temperature. Using permeability coefficient 
p w = 100^m/s and assuming that water enters the complex only via permeation and spreads 
evenly into all water layers, we estimate that permeation causes each water layer to swell 
at the rate of d w ~ 2 x 10~ 3 nm/s, or equivalently, 1.8% of the equilibrium spacing per 
minute (FigEl). Given the disk diameter = 15/xm and the swelling rate d w , the flux of water 
permeating into the disk part of the complex is 

Jd ~ 0.06/um 3 /s, 

which only accounts for a small portion of the total water influx 

Jt ~ l/xm 3 /s 

inferred from the observed myelin growth rate ~ 0.5/im/s (FigEJ- The growth of myelins 
implies that both the bilayer area A and volume V of the complex increase with time. Our 
model, nevertheless, requires that the concentration A/V is maintained at a sufficiently high 
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level so that the disk is always unstable against the formation of myelins. Since Jd <C Jt, 
the effect of water permeation is too weak to cause any significant decrease in A/V and thus 
to compromise the proposed mechanism for myelin formation. By the same token, the flux 
of water permeating into a myelin in Figure [21 is given by 

J m an 2.9 x 1(T 3 L /xm 3 /s, 

where the myelin length L is in jim and the myelin diameter is set to 2.5/im. Now we can 
define a critical myelin length L c using the equation 

Jt = Jd + Jm(L c )- 

When the myelin length L > L c , our mechanism for myelin formation no longer holds be- 
cause the flux due to water permeation, Jd + J m exceeds the total influx Jt and thus the 
concentration A/V must decrease, restoring the stability of the disk morphology. Given Jt, 
Jd and J m calculated above, L c is about 320^m. In Figure El the lengths of all the myelins 
are less than L c , which means our model for myelin formation can be applied to this system 
with self-consistency. 

6 Discussion 

The model presented in Sections 01 and [U offers a simple geometrical explanation for the 
formation of myelins in our experiment: The bilayer repeat spacing increases and therefore 
the inter-bilayer repulsion decreases when a multilamellar disk is transformed into a myelin 
under the constraints of fixed volume, bilayer area, and number of bilayers (see (jllj) and 
(|12(0 . If the lipid concentration is sufficiently high, the decrease in the inter-bilayer repulsion 
energy is greater than the curvature energy of the myelin, and thus the disk is unstable against 
myelin formation. Our model can be thought of as a minimal model for myelin formation, 
from which more sophisticated models can be constructed. For systems with non-negligible 
disk perimeter and myelin end-caps, the proposed geometrical mechanism still holds, but the 
expressions for Dd and D m are no longer as simple as (fTTjl and (|12|). respectively [3"Q. In the 
following we will discuss the limitations as well as the implications of our model. 

In Section 13 we impose two artificial (i.e., non-physical) constraints on the myelin ge- 
ometry in order to illustrate the geometrical mechanism that destabilizes the disk: The 
constituent bilayer cylinders of a myelin are uniformly spaced, and the radius Ri = D m /2 
(equation ("""J). With these two constraints the calculations of the energies (0 and (fTTT^ are 
greatly simplified. Freeing these two constraints would not weaken our proposed mechanism 
for myelin formation in that allowing D m to be non-uniform or Ri to vary can only lower the 
myelin energy further. We will investigate the myelin structure without these two constraints 
in another work|31j. 

In our model we neglect the thermal undulations of bilayers [2*31 124j (Fig[""J). The reasons 
are twofold: First, the persistence length:32l "SB] of DMPC bilayers is much larger than 
any lengths of the disk- myelin complexes observed in our experiment (FigEJ. This means 
the bilayers are stiff, and thus their undulations should be negligible. Secondly, although 
bilayer thermal undulations decrease the bending stiffness k c © 33] 1341 135" and increase 
the inter-bilayer pressure P ([7]) . these effects, however, would only help destabilize the 
disk and therefore would not invalidate our model. Moreover, the threshold pressure Pth 

1 The thermal undulations of bilayers alter the functional form of the hydration pressure Pi lf""ll l28l 
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defined in Section [3] is insensitive to bilayer undulations, as shown below: The decrease in 
the inter-bilayer repulsion energy due to the disk-to-myelin transition can be approximated 
with P th A(D m - D d ) = P th AD d /(2N - 1) (see (JTDJ) and (Q"""J)). B Y definition P th satisfies (see 
©-©, and P1Q 



Using N = 200 and equation @, and assuming D d is about equal to the equilibrium value 
of 6.22nm for DMPC bilayers (Section [""J), the above equation yields the threshold pressure 
P t /j w 0.3atm regardless of bilayer undulations. This estimate of Pth is close to the result, 
0.32atm, obtained in Section In the above calculation we neglect the weak dependence of 
k c on bilayer undulations "*"""! 13*5] . 

Based on the experimental evidence described in Section [21 we postulate that the bi- 
layers can exchange materials through defects to reach quasi-equilibrium (postulate (b) 
in Section |3J. This postulate is important in our model in that it greatly simplifies the 
problem. Therefore the formation of defects in systems like ours deserve more detailed 
investigation [n"J [""3 """5] . 

Our theory says nothing about how the myelin size is determined. This size appears to 
be determined by pre-existing structure in the dry lipid from which the disk and myelins 
grow (Fig|2j). However, our theory does suggest consequences of changes in size: Since the 
bilayer repeat spacing is of the order of the equilibrium spacing, the overall diameter of the 
myelin tube is proportional to the number of bilayers N. We see from equation (|15j) that the 
threshold pressure Pth is proportional to In N/N for large N. Thus an increase in the myelin 
diameter should produce a roughly proportionate decrease in the internal pressure. 

Our analysis of inter-bilayer interaction energy assumes that the bilayers are flat. How- 
ever, in the myelin structure the bilayers are curved. This in principle induces a correction 
to the interaction energy per unit area In the following we show that this correction is 
negligible when the number of bilayers N 3> 1: Equations (["" ) )— lf"" |> . © an d l[14|) imply that 
the net change of the interaction energy due to the spacing increase is approximately 



Using equation (jHJ) the interaction energy density (fTU|) of two adjacent bilayer cylinders of 
radii ~ R in a myelin is about 



where the dimensionless prefactor c is expected to be of order 1. The correction term 
c(D m /R) 2 results from the bilayer curvature 1/R. Terms of order (D m /R) 1 would depend on 
the direction of curvature, and are thus ruled out by symmetry. With ("5J), (JJJJ) and (fT7|) the 
correction to the inter-bilayer interaction energy of a myelin due to nonzero bilayer curvatures 
is given by 




A(iV-l/2) 2 ln(2iV-l) 
Dj-N 3 (N + 1/2) 



(15) 



P x A{D m - D d ) ~ -LDl^N ~ N. 



(16) 




(17) 




N 




71=1 



ln(2iV) 



(18) 
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where the decay length A < D m (JBJ), and R n = Ri + (n — l)D m and A n = 2irR n L are the 
radius and area of the n-th bilayer cylinder of the myelin, respectively. The above estimates 
suggest that the curvature effect (fT8]l is negligible compared to our main effect (|T6|) when 
iV>l. 

In Section we show that the repeat spacing increases when an iV-bilayer disk is trans- 
formed into an iV-bilayer myelin (see and (|12j0 . The spacing would increase further if 
the myelin were converted into concentric spheres or onions 4 also composed of N uniformly 
spaced bilayers. Specifically, the onion spacing ~ (1 + 1/(2N — l))D m > D rn for N 3> 1. 
This suggests that our model might help understand the formation of onions observed in [I] 
and [2]- In our experiment (FigEJ their formation is, however, kinetically blocked because 
they cannot form continuously from a disk or myelin. Furthermore, the relative stability of 
onions is also influenced by the Gaussian bending stiffness |19|. since the transformation from 
a disk or myelin to onions changes the topology 27 . 

Although our model is invented mainly to account for the disk-to-myelin transition shown 
in Fig|2 it also sheds some light on the formation of myelins observed in contact experiments 
(see FigQ) 13 El : In a contact experiment water is brought into contact with concentrated 
surfactant. Immediately after contact, the surfactant molecules self-organize to form many 
multilamellae at the interface [IS]. These multilamellae should have high lipid concentrations, 
and therefore their inter-bilayer interaction is strongly repulsive. As a result, the bilayers 
along the contact interface would curve to form myelinic structures. 

7 Conclusions 

In this paper we propose a geometrical mechanism to account for the formation of myelins: 
In short, if a stack of flat bilayers is internally confined and the inter-bilayer interaction is 
repulsive, geometrical packing alone can lead to myelin formation. We believe this geometric 
mechanism can help explain why myelins form in a variety of experiments El IS]) where 
internal confinement and inter-bilayer repulsion also appear to play important roles. Our 
findings may help develop techniques in controlled growth of myelin or onion structures, 
which have potential applications in encapsulation and drug delivery We emphasize that 
our model, by its very nature, addresses only the equilibrium aspect of myelin formation. 
Many questions such as the myelin size distribution and the myelin growth rate still await 
answers. 
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